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1. Introduction 



The properties of QCD at a finite temperature T and chemical potential n play an 
important role for the physics of the Early Universe, heavy ion collision experiments, 
and neutron stars. Despite the presence of a potentially large scale, say T 3> Aqcdj 
QCD however remains sensitive to non-perturbative infrared physics under these con- 
ditions [|l], 1^ . This means that various effective theories may be derived systematically, 
but the rich phenomena encoded in them (for a recent review, see |0) may have to be 
addressed non-perturbatively. 

Quite generically, effective theories tend to possess extra symmetries, only broken 
by higher dimensional operators. A familiar example is weak interaction induced 
strangeness violation in zero temperature QCD. 

Similar phenomena take place also in high temperature physics. For instance, the 
dimensionally reduced effective field theory describing the thermodynamics of the 
electroweak sector of the Standard Model has extra symmetries: parity (P) and charge 
conjugation (C) are only broken by higher dimensional operators 0. 

In QCD, perturbative interactions break neither P nor C, but the latter can be 
broken by the thermal ensemble, if there is a finite baryon density (yU 7^ 0). Indeed, in 
the dimensionally reduced effective field theory for the thermodynamics of QCD, there 
is again a new C odd operator, of higher order than the leading C even operator . 

The purpose of this paper is to address the same phenomenon in the context of the 
so called Hard Thermal Loop effective theory (for a review and references, see 0), 
the generalisation of the dimensionally reduced theory for real time observables. We 
show that a purely gluonic C odd operator is induced by quark loops. Thus C odd 
observables may be directly sensitive to non-perturbative bosonic dynamics. 

The environment in which we may envisage our results to have significance is mainly 
that of heavy ion collision experiments. In cosmology the baryon chemical potential is 
too small to have any significance, fi/T ~ 10~^. In neutron stars, on the other hand, 
the interesting part of non-perturbative QCD dynamics appears to be related more to 
quark pairing near the Fermi surface than to gluons [Q. 

As far as heavy ion collisions go, precision studies of dimensional reduction show 
that effective theories of the type considered may be quantitatively accurate down to 
T ~ 2Tc IP, |, 1^, [10], [Tl|. At the same time, the infrared dynamics described by the 
effective theories is completely non-perturbative and does not allow for a perturbative 
treatment at any reasonable temperature ||, ||, |12|, |T^. (Apart perhaps 
from observables such as the free energy density where non-perturbative effects happen 
to almost cancel out numerically [0.) Direct four- dimensional (4d) lattice simulations 
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are not available either for real time quantities (for a review, see |TB|). Thus effective 
theories appear presently to be the only way of studying quantitatively some interesting 
non-perturbative processes for phenomenologically relevant temperatures. 

Among the C odd processes one could think of are that various correlation lengths, 
decay rates and oscillation frequencies change, because previously distinct quantum 
number channels couple to each other in the presence of 7^ Ijl^. In particular the 



fluctuations of baryon and energy densities are correlated. One could also consider 
direct C odd observables, such as M = J2tt±{p+ ~ P-)/(P+ + P-) EH Ull) where the 



sum is over the momenta of all charged pions. In addition one may consider contribu- 
tions from C odd amplitudes to C even observables such as dilepton production ^0 . 
We do not here study any of these applications, though. 

The plan of the paper is the following. We fix our basic notation and conventions 
in Sec. ^, review briefly the effective theory in the static limit in Sec. and present our 
derivation of the general non-static case in Sec. H. In Sec. || we argue that the result 
of Sec. ^ is equivalent to another effective description, "classical kinetic theory" , which 
also allows for a simple numerical lattice implementation. We conclude in Sec. |[ 



2. Conventions 

The partition function of QCD at a finite temperature T and chemical potential fi is 

Z = Tre-^(^-^^) = j VA^Vi)Vi)e^Y>[-\sE), (1) 
Se = dr Jd^xCE, (2) 



jO-e = ^Tr F^^Ff,^ + tp['y^D^, - 70/^]^/^, (3) 

where (3 = T~^, H is the Hamiltonian, Q is the quark number operator (three times 
the baryon number operator), boundary conditions over the time direction are periodic 
for Af, and antiperiodic for tp,^, = A^T", Tr T"T^ = 5"V2, = d^- igA^, 
F^u = ii/ 9)[Dfj,i F)y\, and the metric is (+ + ++). We set ^ = 1 in the following. As we 
see from Eq. (|^), the chemical potential corresponds in momentum space to shifting 
fermionic Matsubara frequencies a;^ as a;^ — > 0;^ -t- z/z. 

The observables we are fundamentally interested in are real time correlation functions 
of the type 

C(t„ Xi- T, /i) = Tr e-^(^-^«) \0^ih,x^)02{t2, x^)..] , (4) 
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where Oj are operators in the Heisenberg picture. As usual, such expectation values 
can be found by computing first the corresponding objects in the Euclidean theory of 
Eqs. and performing then an appropriate analytic continuation. 

In view of the analytic continuation, we will at a number of points already write the 

action in Minkowski notation corresponding to the metric g^'^ = diag (H ). The 

continuation will be understood to be made by writing r = it,dr = —idt, uju = —iuj, 
Aq = —iA^. An additional minus sign is inserted in the relation Ce{t = it) = —Cm-, 
such that 

— Se = — J drcfxCE ^ y dtd^xCM = iSu- (5) 
In the continuation scalar products change as 

aE-bE = a^b^ -aM ■ = -df^b^, (6) 

where we use the implicit notation that both indices down implies Euclidean metric. 
We denote the various integrations arising as follows: 

J = Jd^xdt; (7) 
dpo 



1^2^) ' = (Po, p) ; (9) 
TJ2, P, = {un,p), (10) 

where Un are fermionic {uD or bosonic (cu^) Matsubara frequencies; and 

dQy 



Air 



v^ = {l,v'), v^v'' = 0, (11) 



where the integral is over the directions of v^- The corresponding 5-functions are 
assumed normalized such that J^^^: = 1. We also denote 

5_,(p2)^2%o)(27r)5(p2), P'^P-P, (12) 

such that 

/ 54n/(Po,p) = / = |P|. (13) 

All our results will factorise in a form where a phase space integral is left over, which 
can be carried out explicitly. Let 

npiup) = , N±{up) = npiujp - fi) ±nF{up + n). (14) 

6 ~r -L 
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Then, irrespective of the relative magnitudes of T, /i, 



/^iv_M^f(r^ + ^). (15) 

N+{ujp) 1 f dN+iujp) 1 d f ^ , 1^T2 fj^ 



iV_(cjp) /• 1 dN_{up) 1 f d^N_{ujp) /• /i 



cjp Jp Up dujp 2Jp dujp 2(9/i 



3. Effective action at the static limit 

We start the actual computation by recalling the fermion induced C odd operators in 
the static limit. This case corresponds to dimensional reduction [Q. 

In Minkowski notation, charge conjugation can be defined for the bosonic fields as 

A, ---a;, d,^d;, f,,^-f;^. (is) 

Furthermore the Minkowski action Sm should be real. Thus one would expect C odd 
operators to have an odd number of gauge fields. In a non-Abelian case no operator 
with a single power exists, and hence the lowest non-trivial order is cubic. 

Indeed, apart from renormalisation effects, the fermionic contribution (with Nf 
flavours) to the effective action of the bosonic Matsubara zero modes is 

The quadratic and quartic terms here conserve C. The linear term breaks C, but is only 
relevant for an Abelian theory where TrAo = Aq; it has been discussed, e.g., in [pi] . 



The cubic term, the first C breaking operator in QCD, has been discussed in [|]. There 
are no higher order operators involving Aq only P, There are other higher order 
C breaking operators, though, such as 

iC', = -.,N,£,^^{l^O(^)y.A,n, (20) 

but their contribution is suppressed at least by C(p^/(27rT)^), where |p| <gT is the 
dynamical scale within the effective theory. Therefore we ignore them here. 

For future reference, we note that if we define = ~{5 / 5A'^)Sm in the non-Abelian 
case, and j'^ = —{6/6A^)Sm in the Abelian, then, after the analytic continuation 
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Figure 1: The Feynman graph computed in Sec. 



discussed around Eq. (H), Eq. (|^) corresponds to having in the static hmit 

-,2 /T^2 ,,2x „3 
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T2 /i^ 



9' 



(21) 

(22) 



where we have left out the contributions from the last term in Eq. (|T^), going beyond 
the analysis in this paper. 



4. The non-static case 

We now move to the non-static case, to generalise the effective action in Eq. (|l^). To 
determine the leading bosonic C odd operator in QCD, we compute the graph in Fig. |1]. 

4.1. The 3-point function 

The computation of the graph in Fig. |I] is a straightforward exercise. Let us express 
the result as a contribution to the Euclidean action in momentum space. Then, for an 
arbitrary gauge field configuration A^{x), we obtain 



6Se = ta'^d Sq+r+sT^u^{Q,R,S)Tt\a^{Q)A,{R)A^{S)], (23) 

3 JQ,R,S L J 



'JQ,R,S 

where 



F,UPo,p) = ^Tr[f7^(f -(^)7,(f +^)7^]. (25) 
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To evaluate the sum over ul, we use the contour formula 



T /(nvrT + z/i) 



E 



i Res f{z) 

gi/Bz+f^fi _|_ ]_ ■ 



n odd poles at Imz ^ ^ 

In addition we denote 

ul = ujI_q = (p - q)^ = (p + s f. 

To consider only the part odd in C, we replace in the whole expression 

1 



(26) 



(27) 



(28) 



Picking up the poles according to Eq. (^), and using momentum conservation as 
well as the identity np^—ujp) = 1 — np^up) (cf. Eq. (p!^), the expression in Eq. 
becomes 



fiua 



NJtup) 



+ 



+ 



N_{ujp+s) 



-lUJp.p) 



p+s 



[(go + it^pf + ^I-q] [{sq - iujpf + iolj^s] 

F^vajgo -iujp-Q,p) 

[(go - ii^p^qY + ^IWq + i^^p-qY + ^^l+s] 
Ff,ua{-So - iujp+s,p) 



[{so + iuJp+sY + ^p] [{^0 - i^P+sY + ^p-q] 



[LOp 



\^p+s 



-Up) 



(29) 



-^P-q) 



-i^P+S) 



In fact each of these terms gives an identical contribution to Eq. (|23|): changing inte- 
gration variables on the second line such that p — p + q (so that cup-g wp), on the 
third line such that p p — s (so that Up+s ^p)) and in the terms obtained with 
{ujp — > —ujp) such that p — p, we arrive at 



fiua 



N_{up) 



F^j,^(-2CJp,p) 



fiua 



lUp, -pj 



+ (^cychc permutations of yU, u, a; Q, -R, S 



(30) 



where we have denoted ve,ii = {—i,Pi/i^p)- But the cyclic permutations are automati- 
cally reproduced by the trace in Eq. (P3|), so we can replace them by a factor 3. 

We now go to Minkowski space as discussed in the paragraph around Eq. (|^), 
and carry out the small effective coupling expansion. In other words, we look for 
the leading term in the expansion in small Q/up, R/up, S/up, where parametrically 
Q,R, S < max{gT, gfi), while the integration variable gets its contributions from ujp ~ 
max(T, yu) (cf. Eqs. (p!5D-(p!7D). Naively, the leading term in Eq. (^) is of order 
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0{{QS)^^). This however vanishes after cychc permutations, due to momentum con- 
servation (cf. Eq. (p^)). There is no term of order 0{Q~^) due to the symmetries of 
the expression, and then the leading term is of order 0{Q^). In order to find it out, we 
need to know F^^^^ to order 0{Q^), because it is multiphed by a term ~ 0{{QS)~^). 

Taking traces over Dirac gamma matrices in Eq. (|25|) (aheady transformed to 
Minkowski space and with indices raised), we obtain 



4a;| 



+ 2uji 



+ Up -vQ {g'^'S^' - g^^'S^ - g^'^S") -v-S (g^'^Q'' - g^^Q^" - g^'^'Q' 



-Q ■ S iv'^g^"' - ifg""" 



(31) 



where now = /up)^ f ■ f = 0. 

The terms odd in cup in Eq. (pT]) pick up the first and third terms in the expansion 
of the denominators of Eq. (|30|) in Q/up, S/up, while the term even in tup picks up the 
second term. Using in addition momentum conservation to simplify the expression, we 
finally obtain 



X 



/ 



^Q+R+s Tr 



A,{Q)A,{R)A,{S) 



N_{ljp) 



2\2 



+ 



2\2 



ivQ){vR)\{vQ)^ {vQ){vR) {v ■ Rf 
v^'v^R'' + v'^v'^R" - v^'v''Q'' - v^v^'Q'' f Q'^ R^ 



{vQ){vR) 



vS 



v-S V ■ R V ■ Q 



+ 2 



V ■ Q V ■ R 



vS 



(32) 



This result can still be made more transparent, however. We write the gauge field as 



A,{Q) = A,{Q,v) 



v-Q 



A^iQ), i^(g,t;)= r 



v-Q 



(33) 



such that V ■ A = 0. It is then straightforward, employing the momentum conservation 
identity 

(34) 



{vQ)ivR) {vQ){vS) {vS){vR) 
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A^{Q) A,{R) A^{S) 



Figure 2: The Feynman graphs discussed in Sec. 4.2 



to show that only terms involving Afj,{Q,v) survive in Eq. (P^. Furthermore, since 
Afj,{Q,v) is transverse with respect to f^, only the last term in Eq. ( P^ ) gives a con- 
tribution. This result is what would have been obtained if we had relied on the gauge 
choice V ■ A = to begin with. Going furthermore to x-space, we arrive at 



6S 



M 



Jp Up Jx,v ■ d 

1 



27r2 



g'f,Nf Tr[i,>(-— 

J X,V '- ^ U • U 



(35) 



where the last integral was from Eq. ([T7|) . This is the C odd part of Fig. |T| at leading 
order in the small coupling expansion, for Q,R, S < max {gT, gn) <^ max (T, /i). 



4.2. Higher point functions 

In a non-Abelian theory, the effective action in Eq. (^) is not explicitly gauge invariant. 
To obtain a gauge invariant result one has to account also for n-point functions with 
n > 3 at leading order in g. They contribute to the effective action at the same 
order as Eq. (|35|), if we consider non-perturbative gauge field configurations such that 
d ~ gA ~ g'^T. Fortunately such higher point functions can easily be computed, if we 
rely on the gauge choice v ■ A = Q ||2^, |2^, as we now show. 

Let us first recall the idea behind the gauge choice v-A = 0. Once we have carried out 
the frequency sum in Eq. (12^), we are essentially evaluating the propagator of an on- 



shell fermion in a background gauge field (see Fig. |^). Because the fermion is on-shell, 
the result should be gauge choice independent for each p separately, or equivalently, 
for each v^. Thus we can choose a gauge separately for each v'^ such that v-A = 0, 
carry out a simplified analysis in terms of A^, and in the end write the result in a gauge 
invariant form, whereby the gauge choice can be removed. 

In this gauge it is easy to see that the only higher point 1-loop graph which can 
contribute to the term linear in /i (such as in Eq. (|35D) is quartic in A. Indeed, after 
the frequency sum, the result for an n-point function is parametrically of the form 
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(cf. Eq. _ 

r f Njup) Tr[4"(g)f-r""^] 

' ^''^ 

where ^ = 7^^^, < m < ra, and Q denotes collectively the soft external momenta. 
Taking the trace one obtains scalar products of A, P and the Q's. Now if m > n/2, 
there must be at least one product or P ■ A. Both of them vanish: the first because 
the frequency sum puts us on the pole (cf. Eq. (pOD), and the second because of the 
gauge choice. Thus the leading p-integral corresponds to m = n/2, and is 

But this integral can lead to a linear dimensionality, i.e. a term proportional to /i, only 
for n < 4. 

We have thus computed the graph for n = 4, 



6Se oc g^NfT Sq+r+s+t T _ 

JQ,R,S,T JP={uii- 



^^l,p) 



Tr 4(g) (f -Q)4{R) if -Q -$ms) ir +miT)] 

P^{P -QY{P -Q - RY{P + Tf ■ ^ ^ 

Taking again into account that = after the frequency sum (cf. Eq. (|30|)), that 
because of the gauge choice, f anticommutes with 4? and that we need the highest 
power of P possible in the numerator, the relevant part of the trace is 

Tr [mQ)QkR)ms)mn = 2{p ■q){p- t)t, [mmRmsmn. m 

But then the integral remaining is (written now in Minkowski space) 

AQ f A f N4^p) / Tr4(g)4(p)4(g)4(r)] 

dbM ^ g Nf Oq+r+s+t / 2 — / tn I D\ • ^^^^ 

Jq,r,s,t Jp Up Jv V ■ [Q + R) 

This however vanishes, as can be seen using momentum conservation and cyclic per- 
mutations of the trace. We thus find that the correct result for the term proportional 
to II, including all 1-loop corrections, is already given by Eq. (P5[). 
The specific gauge choice made can now be relaxed by writing 

i = -J—v^F^., -^ff'A^ = -^V^—v^F^, (41) 
^ vV ^Vd v-V vV ^ ^ ^ 

where V is the covariant derivative in the adjoint representation. The right hand sides 

of Eq. (^iD are gauge covariant, and taking the trace in Eq. (|35|) makes their product 

gauge invariant. Thus in an arbitrary gauge we get 

SSm = -^,9'^^Nf f Tt{-^v"F^,){-^v('F^'^){-^V'^^v^F,^). (42) 
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This is our final result for the leading bosonic C odd operator, to be added to the usual 
Hard Thermal Loop action, given in ||2^ . 

We have not carried out explicitly the angular integral in Eq. (H^l). Various such 



integrals are discussed in |26|. The structure of these integrals is rather non-trivial and 
it is, for instance, not at all obvious that Eq. (|^) reduces to the third term in Eq. ( |I9D 
in the static limit, without carrying out the integration explicitly. We shall return to 
this issue in the next section, where we find a much simpler way of checking that the 
correct static limit is reproduced. 



5. Classical kinetic theory 

In Sec. ^ we computed the leading bosonic C odd operator in the Hard Thermal Loop 
action, Eq. (^). The form of this operator is however not particularly useful for any 
practical applications. Let us therefore show that the same physics is contained in a 
much simpler description, that of classical kinetic theory (for a review, see p8|). 



For notational simplicity we shall work within an Abelian theory for the moment, and 



return to the non-Abelian case only in Sec. 5.5 



5.1. Setup 

We follow here the pedagogic presentation in Ref. |^ (based on the original work 



m 



2^, ISO, 1^). Let us start with classical electrodynamics, and define 



Then the coUisionless Boltzmann equation for hard particles in a gauge field background 
becomes 

%^-P"fl^+«^/|^)-0. (44) 



dt V5x° " dpf^J 

Let us note that we may in general assume that 

f{x,p) =S+{p^)f{x,p), (45) 

where 5+(p^) is defined as in Eq. (0), since this form is conserved by Eq. (|4^) , due to 

Ki^/^5(p') = 2py^F^f,5\p^) = 0. (46) 

The derivate of 0{pq) included in 5+(p^) can be safely ignored as well, since it would 
contribute only at the point po = p = 0, and has no effect after integration over p. 
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We formally solve Eq. (0) in powers of gF^,,: / = /o + /i + /2 + ••• • This leads to 
the recursion relation 

P ■ dU,{x,p) = -gp-F^,{xf-^^p^. (47) 

OPf3 

The zeroth order gives 

p-dUx,p) = Q. (48) 

We take as a solution a space-time independent function depending, in view of Eq. (^), 
non-trivially only on pq, parametrized by T, /x, and applying separately to all particle 
species i: 

fo^ = hip') fi'\po; T, /i.) = 5+(/) n^ipo - /i.), (49) 

where nF(po) is in Eq. (|I4D . Furthermore, antiparticles are always assumed to come with 
the opposite signs of g and n than particles. Thus, a single Dirac fermion contributes 
two degrees of freedom with +g, two with —g, — /i. 

In addition to these equations, we need the definition of the current induced by the 
hard particles: 

f{x) = -Y.9^ I P'f'^ix.p). (50) 
The equations of motion are 

"^^^ = " I i^M^^^^ Ja/m = duF^'' = f. (51) 

The expression for in terms of the background gauge field thus implies a non- 
local effective action 5*^/ = + 6Sm for the gauge fields only, where 6Sm is to be 
determined from 

^ 5SM = -f. (52) 



5.2. Linear term 

Let us now work out explicit expressions. We start by considering / = /o- Summing 
over Nf flavours, each with two degrees of freedom with +g, +/i and two with —g, —fi, 
we obtain from Eqs. (|49|) , (|50|) , 

f = -2gNf l^5+ip')p^N4po) = ~5>^' gNfi^(T' + ^), (53) 
where we used Eqs. (p^, (p!4D, (p^Sf). Eq. (|5^) is then trivially solved. The result. 



6S 



M 



I (54) 

J X 



reproduces the first term in Eq. ([19|) after analytic continuation. 
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5.3. Quadratic term 



For the next term the considerations are well-known, but for completeness we briefly 
present them here, too. We need 



Inserting into Eq. (|50|) , 

As described in [^, this corresponds to the usual Hard Thermal Loop action |^ 

SSm = -g'N, I ^^±^ I (57) 
jp Up Jx,v [V ■ a)'' 

To check this result at the static limit, it is actually convenient to start from the 
form in Eq. (^) rather than that in Eq. (|57D. We shall discuss this is detail for the 
non-Abehan case in Sec. |5]^, and following that argument, it is easy to see that the 
second term in Eq. ( p2D immediately follows from Eq. (0). 

5.4. Cubic term 

We then proceed to the next order in gF^^. We are not aware of previous analyses 
going beyond Eq. (|57|) in this way. 

To proceed, it is useful to Fourier transform with respect to x: 

f{x,p) = [ e^« V(Q,p), f{Q,p) = I e-^« V(a:,p). (58) 

Solving for /2 from Eq. ( ^TD and then inserting into Eq. (|SD|), we obtain 



R,S Jp ip- Q opp 



ip ■ S dp. 
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T.9f f I 5Q~R~s^-i^^(^)p''F^,{R)-^p'Fs,{S) 
, Jr,s Jp op^ op3 \tp ■ Q J ip ■ b 



dp-y dpp \ip ■ Qj ip ■ S 



_dpii \p-Q 

where we carried out two partial integrations and wrote (cf. Eq. (|4l|) ) 

, d rp'' 
dpf^ ^p ■ R 



(59) 



p'^F^(,{R)=t{p-R)Af,{R)=t{p-Rf^{^)A,{R). (60) 
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We can now take the partial derivatives d / dp^,d / dpp, using 

dA^{Q) Q, - 



P-Q 



(61) 



following from Eq. (|HD| ) (or Eq. (P^D). By a change of summation and integration 
variables we can also symmetrise the expression with respect to (i? ^ S", /3 7), since 
Af3, commute. Using in addition that 



(62) 



we arrive at 



g'Nf 



R,S JP 



5q-r-s hiP^)N4po) 



p ■ R p ■ S 



d 



dp"" \p ■ Q 

Let us now show that the same expression is obtained by computing 



Af,iR)A,{S). (63) 



6 



SAf,{x) 



6S 



M 



6 



-6S 



Ml 



(64) 



where we used 5A^{Q')/5A^{x) = g^^e and we expect (cf. Eq. (|35|)) 



SS 



M 



Q,R,S JP 



OQ+R+S s+{p') iv_(po)i.(i?)i"(5)-^i,(g). 

P ■ V 



(65) 



Indeed, taking derivatives using the right-most side of Eq. (^) and writing the result 
again in terms of Aj3,A^, we immediately recover Eq. (|63D. Thus the C odd operator 
in Eq. (|42D is reproduced by classical kinetic theory, at least in the Abelian case. 

It is interesting to note that to reproduce the action in Eq. (^), we had to carry 
out two partial integrations in Eq. (|59|). In contrast, the static limit, the last term 
in Eq. (|22|), is much more easily obtained from the first row of Eq. (p^. We again 
postpone the discussion on this to the non- Abelian case in the next section. 



5.5. Kinetic equations for the non-Abelian case 

To complete the discussion in the previous sections, we review here briefly the form of 
the kinetic equations in the non-Abelian case, and show that they reproduce the static 
limit discussed in Sec. ^. Because of a considerable proliferation of formulae, without 
any additional physical insight involved as far as we can judge, we do not here present 
a complete non-static analysis, though, as we did for the Abehan case. 
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The simplest way to display the non-Abelian kinetic equations is to follow Eq. (^), 
but replace / by an N^. x N^. matrix. The coUisionless QCD Boltzmann equation (for 
finer details of terminology, see |2^, |2B[) for each single fundamentally charged fermionic 
degree of freedom, and the corresponding gauge current induced, can be written as 



P 



■A/]+fKF,.,^}^o, 

= -g ( p^Tr \T-f . 



(66) 
(67) 



To express this in a more familiar way, we may write the matrix / in terms of a singlet 
distribution function / and an adjoint (or octet) distribution function 



/(x,p) = ^/(x,p) + 2T'^r(x,p). 



(68) 



Projecting then Eq. ( |66D with Tr[...] and Tr[T"...], the equations obtain their usual 
forms (with our sign conventions) 



p-df + gp^P 



dpu 



31 



J:9^lp,r^'^■ 



2Nr 



dp,. 



(69) 

(70) 
(71) 



When computing the current, each quark flavour now comes with colours, in addi- 
tion to two spin degrees of freedom, both for particles and for anti-particles. 

The equations can again be solved iteratively in gF^j^^y-. / = /o + /i + /2 + ••• • At the 
zeroth order, 



Jo 



a{i) 



5+{p^) npipo - fii 
0. 



Iterating this, we obtain at the first and second order. 



{p ■ vrfi^'^ 

P ■ dfT 

{p ■ vrif"^ 



0, 



-^p^F;o5+(p2)<(po-/x.), 



8f 

-av^F" 



a(i) 



dff^ 
dp, ■ 



(72) 
(73) 



(74) 
(75) 

(76) 

(77) 
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In order to argue that these equations contain the same physics as Eq. P^), we 
shall complement the full proof for the Abelian case in Sec. p.4| , by showing that the 
full non- Abelian static limit of Eq. (^) is also reproduced. To demonstrate the latter 
point, we start from the first order in gF^^. In Eq. (75) we can write 



p^F^o = {p ■ VAoT - doip^A^) = {p ■ vr'Al (78) 

so that 

fl^'^ = -^A^oS^ip')n'APo-^^^). (79) 
Inserting into Eq. ([7l|) , we obtain 

r, = -E9^ [ Pjf'^ = 9'NjA^, I p,54p')Kipo) 

= 5,og^NfAl I at; (cup). (80) 

After use of Eq. (|l^) , this agrees with Eq. (^) . 
The equation for /2 , on the other hand, becomes 

{P-Wf2^'^ = ^f"'p'Fl,Al6^{p')n'^{p,- ^.^) 

- ^'^ -d^^p ■ VA^i'Al 64p') n^Po - (81) 



4Ar, 



Here 



so that 



(f'^'ip ■ VA^fAl = ^(p ■ vr\d''''A''X), (82) 



/2^^^ = J-rf"''^4AS5+(p2)<(po-/i.)- (83) 



Carrying out the sum over i in Eq. (|7T|), we finally get 

3l = --/Nfd^'^AlAlj^p,5^{p')N'L{po) 

= -5,,-^g^Nfd'^'^AlAlj^N'L{ujp). (84) 

Using Eq. ([T7|), this agrees with Eq. ([2T|). 
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5.6. Possible simplifications of the classical description? 



We have argued that classical kinetic theory should reproduce the results of the Hard 
Thermal Loop effective theory even for C odd observables, thus correctly representing 
the infrared physics of the system. It has however some merits over the Hard Thermal 
Loop action: it is local, directly Minkowskian and, involving only the solution of classi- 
cal equations of motion, can be implemented numerically without the need for analytic 
continuation. Moreover, judging from non-perturbative studies of dimensional reduc- 
tion, such numerical results should be at least qualitatively reliable for T ~ (2Tc, oo). 

In full generality, the classical distribution functions f{x,p) depend on (4-1-4) coordi- 
nates. If one wants to proceed towards numerical implementation, it is best to reduce 
the dimensionality of the phase space. As we have discussed, f{x,p) = 6^{p'^)f{x,p), 
and the dependence can thus be reduced to, say, the spatial components p. This (4-1-3) 
dimensional problem can already be managed on the lattice [0. However, for the 
usual C even case one can carry out a further simplification, by writing in Eq. ( [TSD 

ff^ = 5+(p') n'APo - /i.) W,yx, v), (85) 



and performing the sum over i and the integral over po iii Eq. (|7lD , reducing thus 
the dependence only to angular variables and a total of (4-1-2) dimensions. Further 



simplifications in the angular variables (like an expansion in spherical harmonics 
0) may also be possible. 

It is then natural to ask whether a similar simplification could be made in the 
presence of /2 . In order to have a proper statistical weighting of the initial conditions 
for the time evolution, one should also work out a Hamiltonian formulation in terms of 
the gauge fields A^, = Fq- and W^, W2, an issue which is not altogether trivial 



3q] . We have not carried out these constructions, but are not aware of any a priori 



fundamental problems in doing so. 



6. Conclusions 

In this paper we have addressed the computation of real time observables in QCD at 
finite temperatures and densities. 

Most of the physical observables one can think of are not computable in perturbation 
theory, due to infrared problems. Real time observables are not well suited for direct 
4d lattice simulations, either. However, it is possible to use perturbation theory to 
construct an effective description of the infrared dynamics. If simple enough, this 
could then be studied by non-perturbative means. 
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We have here completed the bosonic sector of the Hard Thermal Loop effective theory 
for real time observables, by computing the leading C odd operator induced by a finite 
density. This operator is of a relatively simple form (Eq. (^2])) and, as we have argued, 
is reproduced by classical kinetic theory, whose equations of motion are local (Sec. |575| ). 

The classical kinetic theory can then be studied on the lattice, as has been demon- 
strated in the case of the sphaleron rate in the electroweak theory |3^, Q , as well as in 
the case of the defect formation rate in scalar electrodynamics |3^. Strictly speaking 
there are still problems in finding a formulation which has a well-defined continuum 
limit , but these may not be important for practical applications, in which one may 
hope to find an intermediate scaling window allowing to extract physical results. 

When the chemical potential is turned on, the further problem arises that the "Hamil- 
tonian" determining the thermal distribution of initial conditions is complex, resulting 
in a "sign problem": standard Monte Carlo methods are not efficient in sampling the 
configuration space, when one is close to the thermodynamic (infinite volume) limit. 
Numerical tests in the static case have shown, though, that in practice one can reach 
large enough volumes before the sign problem starts to reduce the signal-to-noise ratio 
in any significant way ^ . 



Thus, it seems that a numerical determination of many non-perturbative real time 
quark-gluon plasma observables may become feasible, using classical kinetic theory. 
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